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Abstract 

In this work we study the Lie group analysis of the equation ut+g (u) u^+Xh (u) = 
which is one of the generahzation of the classical invicid Burgers equations (A = 0). 
Seven inequivalent classes of this generalized equation were classified and many exact 
and transformed solutions were obtained for each class. 



1 Introduction 

Damping can be regarded as an effect that tends to stabilize or reduce other effects in a given 
system such as friction in a mechanical system. Damping devices are used in almost all the 
systems that we know, thus the importance of having an equation that models appropriate 
damping systems cannot be overstated. In this paper, we shall use Lie symmetry to perform 
analysis on a generalized invicid Burgers' equation with damping in the form 

Ut + g (u) + Xh (u) = (1) 

where A 7^ is a constant and g {u) and h{u) are smooth functions of u = u{x,t) in the 
domain of definition. Numerous applicable cases of equation ([1]) are available in the literature, 
for instance, researchers have used it to model the Gunn effect in semi conductors, rotating 
thin liquid films, chloride concentration in the kidney and flow of petroleum in underground 
reservioirs, see [7] for all the references. If h (u) = O {u'^), a!>0, 0<m<<1, Murray (see 
[0]) shows that a finite initial disturbance zero outside a finite range in x decays under certain 
conditions and for each condition, the asymptotic speed of propagation of the discontinuity 
was given together with its role in the decay process. Also for X = 1, g = u and h = u{u — 1), 
equation ([1]) becomes the limiting case of the Burgers'-Fishers equation 

Ut + UUx + U{U-l) = -Uxx 
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which has been used to model many physical phenomenon in Mathematical Biology and Ge- 
netics. 

Symmetry analysis based on local transformation group is one of the most powerful and 
prolific methods used for solving nonlinear partial differential equations (PDEs). Its ap- 
plication to the study of PDEs was laid down by Sophus Lie, a Norwegian Mathematician, 
in the later half of the nineteenth century. A symmetry of a system of differential equa- 
tions is a transformation that maps any solution to another solution of the system. Such 
transformations are groups that depend on continuous parameters and consist of either point 
transformations, acting on the systems' space of independent and dependent variables, or, 
more generally, contact transformations, acting on the space of independent and dependent 
variables as well as on all first derivatives of the dependent variables. Lie showed that the al- 
gebra of all vector fields (infinitesimal generators) that leave a given system of PDEs invariant 
could be found by solving over-determined (or under-determined) auxiliary system of linear 
homogeneous PDEs, known as the determining equations of the group. The calculations in- 
volved in obtaining these determining equations and their possible solutions are cumbersome, 
hence researchers have written codes for implementing Lie's algorithms in various packages 
for symbolic computations like Mathematica, Maple, etc. Though these softwares help in 
circumventing the tedious tasks encountered in obtaining similarity solutions, they fall short 
of providing any tangible results when the underlying PDE has an arbitrary function like 
that in ([1]). The key feature of a Lie group, which makes it very useful, is the parametric 
representation of smooth functions on a continuous open interval in e; the group parameter. 
This ensures that the mapping is different iable and invertible and that the mapping functions 
can be expanded in a Taylor series about any value of e. For detail ramifications on the 
subject, see the monographs [2]- [5]. 

This work is organized as follows. In section two, we carry out the group analysis of equation 
([T]) and obtain seven inequivalent classes where both g (u) and h (u) are restricted to order of 
powers of u. In section three, we present theorems on transformed solutions for these classes 
and show how they can be use to generate non-trivial many-parameters invariant solutions 
of equation ([1]) from the trivial solution u{x,t) = 0. Section four deals with invariants 
solutions. while conclusion is presented in section five. 



2 Group Classification 

The generalized invicid Burgers' equation with damping ([1]) is said to be invariant is invariant 
under the first order prolonged symmetry operator X^^'^ = + + rj^^ + V^^;- + 
if and only if 

X^^^[ut + giu)u^, + Xh{u)]\i^=0 (2) 

where 

= Vx + iVu - L) Ux - r^Ut - CuUl - TuU^Ut, 
= Vt-^tUx + {Vu-^t)Ut-^uUxUt-T^u'^^. 

Expanding equation ([2]) leads to under-estimated system of partial differential equations in 
£^{x,t,u), T{x,t,u), and (f){x,t,u) which are the coefficients of symmetry generator to be 
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determined. Due to this predicament, group classifications of equation ([2]) is very difficult, 
if not impossible, to achieve. To circumvent this, we shall assume — '''u = ^ which , when 
applied in the expansion of equation ([2]), results in the following equations 

Xh' (u) (/) + Xh (u) Tt + Xg (u) h (u) — Xh (u) (p^ + (p^ + g (w) (f)^ = 0, (3) 
g' {u)<i)-i^-g iu) i, + g {u) Tt + g"" {u) = 0. (4) 

Differentiating equation (jl]) twice and re-arranging, we get 

^ {g\u)A{x,t)+g{u)B{x,t) + C{x,t)) (5) 



g' {u) 

subject to the consistency conditions 

= -A{x,t), (6) 

C-T, = B{x,t), (7) 

it = C{x,t), (8) 

where A{x,t), B{x,t), and C {x,t) are arbitrary functions of x and t. Using equations 
dSJ — dHD into ([3]) leave us with a very complex equation which is absolutely difficult to solve 
if g and h are left as arbitrary functions of u {x, t). But much can be achieved if both g and 
h are considered as functions of order of u which lead us to the following results. 

Theorem 1 Let g (u) = , h (u) = where k and m are real numbers such that k ^ m 
and both not equal to zero. Then we have the following exclusive (inequivalent) cases for k. 
Case 1. For k ^ ± (m — 1), k i^^^^, k ^ ^^yi and m ^ 1, equation ([I]) admits three 
dimensional symmetry algebra spanned by the closed vector fields 

Ml = a,, M2 = dt, M3 = ^'"^^^ xd^ + ^—j^idt + ^udu. 

Case 2. For k = m — 1 {m ^ 1) , equation admits eight parameters group of projective 
transformations spanned by the closed vector fields Mi,M2, 



tdx + ( - XtU ] du 

m — 1 



^-\{m-l)x 



M3 = td^ + Xim- 1) t^dt + ( 2Xtu \ du, M4 

\m — 1 / 

(—1 7;*" \ 7; 
Y7 T^dt + 7du e^("^-l)^ Me = -tdt + -9,, 
A [m — 1) m — 1 J m — 1 

Mr = ( TT— + —^du^ g-A(m-l)x^ 

\ A [m — 1) m — 1 J 
Ms = + A(m- - A'(m- l)t7z™a,) e^(™-^)". 

Case 3. For k = 1 — m {m ^ 1) , equation ([I]) admits eight parameters group of projective 
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transformations spanned by the vector fields Mi,M2, 



Ma 



+ 



1 



(1 - m) 



-t^ \d^-\xt + 



A(l - m) 



1 — m 



^« + { - m)t' - X ] u + X' {1 ~ my fu" 



3X{1 ~ m) xt 
1 



A' (1 - m) 



-t'']d^+[^Xil-m)t^-x]dt 



1 — m 



u'-"" + ( 3A (1 - m) X - -A' (1 - mf t^ ] u 



Ms = 


|^A(1 


Me = 




Ms = 


— 


Case 4- 


For k 



t^ -xt]da;- t^dt + 



1 — m 



tu+ [ ^X{1 - m)t^ - x ] u''' 



du, 



u 



1 — m 



du,Mj = X 



A(l -m) 



1 — m 



{u + X{l-m) tu"')du, 



m—l 



m 7^ 1), equation admits eight parameters group of projective 



transformations spanned by the vector field Mi,M2, 



M, 



xt — 



A(m - 1) 



x^\d,+ \e + 



A' (m - 1)^ 
16 



d, 



m — l 



X^ (m—l) n m+l 

^ ■ —X-^U 2 



t + -X{m — 1) x"^ ] u + tu 



3 — m 
2 



du, 



Ma 



t - -A (m - 1) ) (9^ - 



A^(m-lf , 



m — l 



-A [m—l) X u X[m — l)xu + u 



3 — m 
2 



<9„ 



Ms 
Mfi 



A 



x^a^ - xt - - (m - 1) x'^ 9t + 



A 



xdx + [ t + ^{m — 1) x'^ ] dt — Xxu 2 du, 
X 



m—l 

m + l 



t + -A (m — 1) x^ ) M 



m + l 

2 — XW 



9„ 



M7 = xOi; H — (m — 1) X c/j H \u — X[m—l)xu 2 ) Mg = — xa^; H u 2 a„. 

2 m — IV / m — l 



l—m 



Case 5. For k 

transformations spanned by the closed vector field 



m^l), equation ([I]) admits three parameters group of projective 



Ml, M2, M3 = \xdx + tdt + —^udu- 
1 l — m 
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Case 6. For k = {m ^ 1), equation ([2]) admits three parameters group of projective 
transformations spanned by the closed vector field 



Ml, M2, M3 = -2xd., - 3tdt 



m — 1 



udu- 



Case 7. For m = 1 (A; 7^ 0), equation ([I]) admits eight parameters group of projective 
transformations spanned by the closed vector field Mi,M2, 



Ms 
Ms 
M7 



Xkx'^d^ -xdt + - + 2Xkxu) du, M4 



-xdt + - + Xkxu) du 



1 

k \ \ 



u 



+ u'-'^du ) e-^^\ Me = xd^ + -d^ 



k \ \ 



dt + udu ] e^"', Ms = Xk{ -xd^ + —dt + Xxu^-^du ] e 



-Xkt 



Proof. Using g = u^, h = and equation ([5]) in ([3]) and simplifying we get 



- (m - k - 1) A + 
k 



-{m-l)B + Tt 



k 

+ y{m + k-l) Cu"^-^ + f (A + B^) ti'^+i + j{Bt + a) u + jCtU 



u 



,1-fc 



k 



k 



k 



0. 



(9) 



Case 1. Since all the coefficient in Q) are independent of u, we have | (m — k — 1) A+r^; = 
which implies ^4 = using equation IQ, thus = and = 0. Also C = implies = 
and Bt = and so B = a^, a constant. Therefore, the required infinitesimal transformations 
are 

k — m + 1 1 — m 1 

C = ; ce^x + ai, r = — ; — a^t + a2, 4> = ycxsu (10) 



k - ■ - 
which proves the result. 

Case 2. Under this case, equation Qj becomes 



k 



^ (A, + A (m - 1) T,) u^^~' + [A (m - 1) (5 + r*) + A + B,] u'^ 



m — 1 



which holds if and only if 



m — 1 

+ (2A {m-l)C + Bt + C^)u+ -^—CtU^-"^ = 0, 

m — 1 m — 1 



C{x,t) = C(x), 
A^. - A(m- 1) A = 0, 
A (m - 1) e,. + A + 5, = 0, 
2A(m- + + = 0. 



(11) 
(12) 

(13) 
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From dHD, (USD and (PS]), we have 

B = -{C' (x) + 2X{m-l)C)t + a (x) , 
A = ^ [C" (x) + X{m-l)C' (x)] - [A (m - 1) b' {x) + a' (x)] t + d{x), 

where a (x) , b (x) and d (x) are abitrary functions. Consistency criterion for A in equations 
(fTTl) results into 



C" (x) - A" (m - 1)' C (x) = 0, 
A^ (m — 1)^ b' (x) — A (m — 1) 6" (x) + X{m — 1) a (x) — a" (x) = 0, 
d' (x) — A (m — 1) (i (x) = 0. 



(14) 



Using (1141) with ([6]) and ([7]), and solving the resulting equations that follow, we get the 
infinitesimal transformations to be 



-X(m—l)x 



a5 



A (m — 1) ait^ 
1 



A(m - 1) 



m — 1 



A(m- 1) 

(ag - A^ (m - 1)^ art) ^x^e^^™"^) 



(-2A (m - 1) ait + a^+ (-A (m - 1) + as) e-^(™-i)^) u 



+ (ai + a2e 



— X{m—l)x'^ y2—m 



The Lie algebras are obtained from these infinitesimals. The closure of these algebras is 
shown in the following commutator table where a = A (m — 1) and V = aM^ — 2Mi. 





Ml 


M2 


M3 


M4 


Ms 


Me 


M7 


Ms 


Ml 











-aM4 







-aMy 


aMs 


M2 








Ml - 2aMfs 


-aMj 





-M2 







M3 





-Ml + 2aM6 








Ma 


M3 


M4 





M4 


aM4 


aMy 








^ + M6 
a " 


M4 





aM'i 


Ms 


-aM^ 





_M» 
a 


-^-Me 





-Ms 


M2 

a 





Me 





M2 


-M3 


-M4 


Ms 











Mr 


aMy 





-M4 





M2 

a 








V 


Ms 


-aMg 


a^Ms 





-aMs 








-V 






Case 7. We have from ([9]), 

A = A{t), 

A! [t) - 2XkA {t) + B^ = 0, (15) 
Xkrt + Bt + a = 0, (16) 
XkC + Ct = 0. (17) 
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Equations (jH]), ffT51) and ffTBl) respectively imply 



r = -y4 (t) X + p (t) , 
B = {2XkA{t)-A'{t))x + q{t), 

C = ^{A"{t)-XkA'{t))x^ -{Xkp'{t) + q'{t))x + r{t), 
where pit), q (t) and r (t) are abitrary functions. Hence equation (|T7I) gives 



A (f) = tti + 026^^* + ae 



-Xkt 



Equations 



Xk ip" (t) + Xkp' (t)) + q" (t) + Xkq' it) 

T8l) and their stability criteria lead to 




a 



Pit) 

lit) 
Bix,t) 

Cix,t) 



'Xkt 



a 



5 Xkt 



aQ + aye — — e 
Xk 

{2Xkai + Xka2e^^*^ x + 04 + a^e 



Xkt 



-Xkt 



ase + iXk) aje'^'^^x 
Thus, the coefficients of the symmetry generator are 

«3 ^-Xkt 



^ = Xkaix^ + (^4 - Xkaje ^''*) 



T = - (ai + 0^26^''*) X - -^e"""- + aQ + aje 



^ Xkt 

Xk 



+ "8, 
Xkt 



1 
k 



{ai + 026^''*) u''+^ + {Xk {2ai + ase^''*) x + 04 + a^e^''^) 
+ (03 + iXk)^ ajx) 



u 



Hi 



The Lie algebras are obtained from the infinitesimals above and their closure is shown in the 
following commutator table, = M2 — 2XkMQ. 





Ml 


M2 


Ma 


M4 


M5 


Me 


M7 


Ms 


Ml 








-iV 


XkMr 





Ml 





A'PMs 


M2 











XkM^ 


-A/^Ms 





XkMj 


-AfcMg 


M3 


iV 











-Ms 


-Ms 


-M4 





M4 


-XkMr 


-XkMi 







JM2 ]\/f 
-TT - Me 


-M4 





-AA;M3 


Ms 





XkM^ 


Ms 


fi + Me 





Ms 


Afc 





Me 


-Ml 





Ms 


M4 


-Ms 











M7 





-XkMj 


M4 





Ml 
Afc 








2M2 - AfcMe 


Ms 


-A'PMs 


XkMs 





AitMs 








-2M2 + A/cMe 






Cases 5 and 6 can be obtained by substituting respective values of k in (fTO!) while the re- 
maining cases can be proved in a similar way as did in case two. ■ 
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3 Group Transformations 

Each of the infinitesimal generators Mj in theorem [1] above can be used to obtain correspond- 
ing one-parameter (e) Lie group of point transformations Gt through exponentiation. To 
achieve this, we solved the following system of differential equations 



dxk (g) 

de 
du{e) 

de 



U {e),u{e)) , 
(t){xk {e),u{e)) , 



subject to the initial conditions 

for all Xi : i = 1,2,. ..,8 and obtained the following (in addition to Gi : {x,t,u) 
(x + £, t, u) and G2 '■ (x, t, u) 1 — > (a;, t + e, u)): 
Case 1. 

fc — m-f-1 



G. : 



Case 2 (a 



Ge : 
Ga : 



Gt. 
Ga : 



Case 3 (a 

Gq : 



Case 4 ( a 

Ge: 
G7: 

Gs : (x, t, u 



A(m- 1)). 



(^x - In ^1 - aet, ^z^, '"Ye (1 - aet) + (1 - aetf u"''^^ . 



ax _\_^,m — l 

'x~ f 5 



X, e~'^t, e™-i?/ 



— > (in v^^FTe^, t, e^^' '"v^-e + e"^n) , 
— > In v^^oFTF^, 
A(l -m)). 



1— m 



(^(x + f t^) e= - f t^, t, i-'^e^ (m^-™ + 2at) - 2at 
(^[x + (1 - e^)] e^ eH, ^"^a (1 - e^) t + li^-^ 



A(m-l) 
2 



e^x, [t + f (e*^ - 1) x^] e^, a (e^ - 1) x + 2 



e=x,t + f (e^^ - l)x2. 



a (e^^ — 1) X + 2 ' ) e 



2 

1 — m 



e X, t, I — £ + M 2 



1 — m \ 1 — m 



Case 5. 



Case 6. 



G3: 
Case 7 

G3 



Gj 
G. 



3 3 

e2^x, e^t, el-"" u 



e ^^x, e ^^t, e™-i^M 



x ^ _ ln(l-AEfca:) u 

1-Aefcx' Afc ' fe/i_Aefcx-ex«fe 



Afc 



ln( — e+e~ 



Afc 



ln(A£fc+e^'=') fc/ A2sfc2a;e-Afct ^ 



1+Aefce-^fe' ' Afc ' V 1+Asfce 



From all the above calculated one parameter Lie groups of transformations, the following 
theorems on transformed solutions hold true. 

Theorem 2 If u = ip{x, t) is an invariant solution of ([I]) obtained through the generators 
Ml and M2, then so are the following functions: 
Gi.ip {x, t) = (p {x — e,t) , 
G2.V2 {x,t) = (p{x,t - e). 

Theorem 3 If u = f{x,t) is an invariant solution of ([I]) obtained through the remaining 
generator of case one, then so is the following function: 

G3.(/? (x, rj = e fey^le ^ x, e k 

Theorem 4 Ifu = ip{x, t) is an invariant solution of ([I]) obtained through the generators of 
case two, then so are the following functions: 



G-i.ip{x,t) 
G^.Lp {x, t) 



1/ [ip(x-ln ^/T3^,^ 



-£{l-aet) 



{1-aety 



1/ {1 + aete-^'=) [ip {\n ^aet+~e^ , t)]"" ^ - ee 



f,ax \ "I 1 — ?Tl 



G,.^{x,t)= ^-^[p[x,t-^-^)] 

G^.ip (x, t) = e~"^ip (x, e~^t), 

Gj.ip (x, t) = e"^^ ""'^/-e + e^'^tp (in ^-e + t) , 



l—m 



a^et 



Theorem 5Ifu = {p{x,t) is an invariant solution o/Q obtained through the generators of 
case three, then so are the following functions: 



■,.ip{x,t) = ' \J[ip{et + x,t)] 



l—m 



Gs.ip{x,t)= ^-'\/[v^((x + f^2(l-e^))e^^e-)]'~' 



at (1 - 



Theorem 6 If u = ip{x, t) is an invariant solution ([I]) obtained through the generators of 
case four, then so are the following functions: 



Ge.ip (x, t) = [if (xe^ {t + f (e^ - 1)) e^)] 
Gr.ip {x, t) = (xe% t + f (e^^ - 1))] 

r 1 — m 

G8.^{x,t) = [^{xe-',t)] 2 +£ 



1 — m 
2 _ 



l — m 
~~2~ 



- ax (e^ — 
ax (e^^ — 1) 



2 

l — m 



2 

l — m 



2 

l — m 



Theorem 7 If u = ip{x,t) is an invariant solution of obtained through the respective 
generator of case five and six, then so are the following respective function: 

Gs.Lp {x, t) = e~ (^e^'^x, e'^t^ . 
G-i.(p {x, t) = e~'^'^Lp (e"^^x, e"^'^t) . 

Theorem 8Ifu = {p{x,t) is an invariant solution of obtained through the generators of 
case seven, then so are the following functions: 



G3.(f{x,t) 

G4.ip {x, t) 

G5.(p{x,t) 



{l—\ekx) 



kj i+ggAfct ^Xkx+ In ^^/Aefcx+e-^*^* j]"*^ 



[vii^rx^t-ln^VT^X^)], 



(x - f-.e-^'^^ 



t)] -ee-^'^'. 



GQ.(p{x,t) = e kif (e^x.t), 

^ 09 ( X, — In £ + e 



GT.(p{x,t) 

Gs-^p {x,t) 



-Xkt 



l+Xeke-^''^ 



, In V Xek + e 



Xkt 



X^k'^exe-^''^ 



Example 9 This example highlights the importance of the above theorems. We start with 
the trivial solution of equation ([I]), that is u{x,t) = 0. By applying the transformations of 
theoremlE on this solution we found that G3, G4, Gq, and G-j all give no new solution but 
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Gs and Gg respectively give u = \J—ee -^'^^ and u = \j ~ '^i^xl^l-\kt as the new non-trivial 
one parameter invariant solutions of ([I]) • We next pick one of these solutions and apply 
the rest of the transformations on it so as to generate more new solutions. Suppose we pick 
u{x,t) = —ee~^^^ which is invariant under G^, the actions of the rest of the Gi are as 
follows: 

Gj:G,^ y^ii 
Gg : Gi : 



Xkt 7 



Lt4 . U-g y (ee2_i)eAfet_Afcei' 



„ „ / {\—\k£-ix){^k'^e\x—\kee-iX^e^ 



„ „ _£4 k (l~\ke3e'^4,x)(^X^k'^eie^4,x-Xkee3e^ix+ej 

Crfi . Lt^ > e k \l ^xkt_^_x(^i^xke3e^ix)(ee3e^4x-kei)+X^k^eie3e'^^ix^ ' 



Hence we see that by applying the transformations of theorem\^to u{x,t) = 0, we obtained 
a non-trivial five-parameters solution 



u[x,t) = e 



(1 — A/cesC^^x) (A^/c^eie^^a; — Xkee^e^'^x + e) 



that keeps equation ^ invariant. If we have chosen u = \J which is invariant 
under Gg, then the actions of the rest of the Gi will be as follows: 

. k / X'^k'^ex 

Lt7 . Lrg y (1-Afc££i)e^'=*+Afe£' 



r~< r~< k £2(Afc££i-l)-A^fc^£a: 

iJ-S • <-f7 Y Afc£+(1-Afc££i)e^'=* ' 



(~i k I £2(Afc££i-l)-A^fc^£3; 

Lt4 . Lt5 \I A/!ce+(l-Afc££i)(l-£2£3)e^'=* ' 



Gs : G4 



J./ (1 — Afc£4) £2(Afc££l — 1)(1 — Afc£4)— A^fc2£zj 



A(l-AA;£4)[fc£-£2£4(l-Afc££i)a;]-A^fc2££^^,2_|.(;^_Afc££i)(l-£2£3)e^'=' ' 



„ „ _£5 fc/ (1 — Afc£4) £2(A/C££l — 1)(1 — Afc£4) — A fc2£e'^5a 

Gq-.Gs ^ e 



A(l— AA;£4)[fc£— £2£4(1— Afc££i)e^5x] — A'^A;2££4e-^^5z2 + (l— Afc££i)(l— £2£3)e^ 



Therefore, 
u (x, t) = 



(1 _ Afe4) [e2 {Xkesi - 1) (1 - AA;6:4) - X^k'^ee^'^x 



A (1 - AA;e:4) [ke - e^e^ (1 - Xkesx) e^sx] - )Kk'^ee^e^''-^x'^ + (1 - Afesi) (1 - £2^3) e^*^*' 
zs t/ie new six-parameters invariant solution of equation 

• Transformations Gi and G2 were not used because they only generate space and time 
translations respectively. 

• The order of applying the transformations is immaterial. 
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• Gi : Gj means the action of Gi on the solution of Gj. 

Theorem 10 // g {u) = h (u) = {k ^ Q), then all the above results hold as follows: 

1. For k ±1, k and k ^ equation ^ admits three parameters group of projective 
transformations spanned by the vector fields of case 1. 

2. For k = 1, equation ([I]) admits eight parameters group of projective transformations 
spanned by the vector fields of case 7. 

3. For k = ^, equation ([I]) admits eight parameters group of projective transformations 
spanned by the vector fields of case 3. 

4- For k = —1, equation ([I]) admits eight parameters group of projective transformations 
spanned by the vector fields of case 4- 

5. For = |, equation ([I]) admits three parameters group of projective transformations 
spanned by the vector fields of case 5. 

6. For k = — J, equation ([I]) admits three parameters group of projective transformations 
spanned by the vector field of case 6. 



4 Reductions and Exact Invariant Solutions 

The reason behind going through all the tedious process of finding symmetry generators 
of any given differential equation is mainly to use them for obtaining symmetry reductions 
and possibly, exact solutions of the underlying differential equation. In this section, we 
shall use the symmetry generators of each of the classified cases in the previous section to 
derive corresponding reduced equations and subsequently, obtain exact invariant solutions of 
equation ([1]) whenever feasible. 

Case 1. The invariant surface condition of the generator M3 = ^~^^^ xdx + ^-j^tdt + ^udu 
is 

k — m + 1 du 1 — m^du 1 
k dx k dt k 
with corresponding characteristics equations 

k dx k dt du 

k- 



k — m + 1 X 1 — m t u 

1 1 1 
Solving these equations give two invariants tp = x''-™-+'^t 1-™ and p = ut 1-'". Since ijj is 

independent of u, then p = F (-0), F arbitrary. So the invariants become 



By chain rule, we have 



1 I 1 



1 ■m- , 

Ut = t—{F-i{jF'), 

1 — m 

1 m — k 

Ux = -, -x^-^+i F 

k — m + 1 
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which changes ([T]) to the reduced equation 

F-^F' + , ^""^ ^^-^F^F' + A (1 - m) = 0. 

This reduced equation is very difficult to solve for arbitrary k and m, but much can be 
obtained when appropriate choices of k and m are made. For example, if m = A; = 2, then 
the above equation becomes 

F-ipF' - F^F' - AF^ = 0. 
After re-arrangement and simple quadrature, we obtain 

- {Xij + A)± J (AV- + Af + 4tp 

F = ^- 

2 

where A is a constant of integration. Hence the required invariant solution of equation ([T]) 
with m = A; = 2 is 



u {x, t) 



{Xxt + A)± J{Xxt + Af + Axt 



2t 

Case 2. The characteristics equations of M3 are 

dx dt im — 1) du 



t A (m - 1) t2 u'^-r'^ _ 2A (m - 1) tu 

1 

having invariants i/j = "^^^""^ and u = (^ x{rn~i)t ~^ j \ Substitution of u and its 

respective derivatives in ([T]) yield the reduced equation tpF' + F = having solution F = ^. 
Thus the required invariant solution is 



A(m— l)x- \ m-i 

' A(m-l)t + t J ■ 

The operators M4 and Mq also produced same result as M3 with the constant of integration A 
replaced with ^^"^^-^ and I'espectively. Reduction with M5 or Mg yield the invariant 

solution (A + A (m — 1) t) 1-™ while u = Ae~^^ is obtained as a solution of ([1]) by using M7. 

Case 3. Two invariants of Mq are ip = t and u = (f + -P) • These lead to the 

reduced equation F' + ^ + A (1 — m) = having solution F = — . Thus the 

solution of ([T]) with /c = 1 — m is 

1 

A + X A (1 — m) \ 

u = 1 

t 2 

Solving the characteristic equations of M-j give the invariants if) = t and u = [{2x + at"^) F — at] 
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with the second one obtained by solving the Bernoulh equation 

du Xu Xatu"^ 



dx a {2x + at^) a {2x + at"^) 

where a = A (1 — m). The reduced equation obtained from these invariants is F' + = 
leading to the solution F = Hence 

^ ^ / 2x + A(l-m) iA-t)t \ ^ ^^g^ 
y 2t — A j 

Reduction with Mg produced same result as that in (I19p with A = 0. Mg also give additional 

1 

solution u= (A — A (1 — m) t) i-™ . 

We next use a linear combination of M-j and Mg to obtain a solution of ([T]). The characteristic 
equations of My + Mg are 

dx dt , du 

— = — = (1 — m) — 
2x t u 

with invariants ijj = u = (tF) i-™ . The reduced equation obtained is 

(F - 2^/'2) F + (F - a) = (20) 
where a = A (m — 1). To solve this innocent looking differential equation , we let 

^ = TW) ■ 



Thus 



F' 



{21/;^ -a)z' + Ai/jz {z - 1) 



z^ 



a-F = '°-^^-'"-'-^>, (22) 

z 

F-2%1)^ = — . 

z 

Substituting (l2T|l and (l22l) into (!20!) lead to the following seperable first order differential 
equation 

dz 2ipdip 
z{z-l){z-2) ~ a -21/;^ 

having solution 

(2 - 1)^ {2ij^ - a) 

Without any lost of generality, we let ^4 = 1. Hence the solution of the reduced equation 
is 

- 2 (a + 1) F + 2^/;^ + a (a + 1) = 
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which is a quadratic equation in F and so 



Therefore 



u 



F= {a + l)±ya + l-2ilj^. 
[A (m - 1) + 1] t ± ^/[Mm^^W+^W^^ 



1 

1 — m 



is the required invariant solution of ([T]) when M7 + Mg is used as reducing generator. 

The resuhs for the remaning cases shall be presented in tabular form, A is a constant 
wherever it appears. 

Case 4 (6 = A(m - 1)). 



Subalgebra 


Invariants 


Reduced Equation 


Solution 


Me 


u= [F+ fx) 


(F - V-) F' 


M = (A + ^x) ^ 


M7 


^ = - t 


F' + 1 = 




Me - iMy 




{2'il;F - + |i?3 _ ^2 _ Q 





Case 5 &; 6. 



Case 


Subalgebra 


Invariants 


Reduced Equation 


5 


M3 




- 3^f) F' + 3F2 + 3A (1 - m) = 


6 


M3 




(3F - 2V/) F' - 2^F + 2A (m - 1) i/jF^ = 



Case 7. 



Subalgebra 


Invariants 


Reduced Equation 


Solution 


M3 


^ (.t2 Afc.x) 


^jF' - F = 




M4 




F' = 


i 


Ms 


ip = t,u = {F ~ XkxY 


F' = 


u= {A - XkxY 


Me 


ip = t,u = (xF)^ 


F' + F2 + AA:F = 


u-{f%-:L)\B-Xke'''^ 


M7 


ip = x,u = [e^^^^F) 


F' = 


u = {Ae-^''^) ^ 


Ms 


i/j = xe^'^^, 
u= {F - XkxY 


F' = 


u= (A - Xkx)^ 



Finally, we present results obtained through the generators Mi = and M2 = dt which 
are known to have travelling waves solution. The invariant of cd^ + dt is x — ct for any 
arbitrary constant c 7^ 0. Since du is missing, we can only proceed by letting u = F (i/j) 
where ip = x — ct, hence Ut = —cF' and = F'. The obtained results are as presented in 
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the following table. 



Case 


Reduced Equation 


Solution 


1 


(i^fe - c) F' + XF"" = 




2 


(^pm-i _c)F' + AF™ = 


^1— nig— cn^ *" X{l—in){x—ct) 


3 


(F™-i - c) 


F' + AF™ = 


= c ± - 2A (1 - m) (x - ct) 


4 


( 2 — c 


F' + AF™ = 


m-1 l±A/l-c(l-m)M-A(x'-ct)] 
U 2 = 5^ 

c 


5 


(F— - c 


F' + AF™ = 


2^/ * 2 ' - 3«i^ = 3 (1 - m) [A - A (z - ct)] 


6 


/ „ m — 1 \ 

(F— - c 


F' + AF™ = 


3m * 3 ' - 2ciii ™ = 2 (1 - m) [A - A (x - ct)] 


7 


{F^ -c)F' + XF = 


y — ck^u'^ Afc(x— ct) 



Remark 1 All the invariant solutions obtained in this section can he subjected to the trans- 
formations of the previous section to generate new solutions. 

5 Conclusion 

In this paper, we have used symmetry analysis to perform classifications and subsequently, 
exhibit many invariant solutions for the damped inviscid Burger's equation ([T]). Though 
various researchers have used this equation to model important physical phenomena, I have 
search the literature but did not see any work that deal with its analysis using Lie method 
as presented in this paper, hence this is the first among the series of work dedicated to the 
analysis of equation ([T]). 
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